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Abstract 

We study the synchronization region of two unidirectionally coupled, in a master-slave configu- 
ration, FitzHugh-Nagumo systems under the influence of external forcing terms. We observe that 
anticipated synchronization is robust to the different types of forcings. We then use the predict- 
prevent control method to suppress unwanted pulses in the master system by using the information 
of the slave output. We find that this method is more efficient than the direct control method based 
on the master. Finally, we observe that a perfect matching between the parameters of the master 
and the slave is not necessary for the control to be efficient. Moreover, this parameter mismatch 
can, in some cases, improve the control. 



PACS numbers: 05.45.Xt,45.80.-|-r,05.40.Ca, 87.19.1r,87.19.1j 



I. INTRODUCTION 



Synchronization of nonlinear systems is a common phenomenon studied in physical, 
chemical and biological systems, among others [T]. One of the most astonishing cases is 
the so-called anticipated synchronization: when two dynamical systems are unidirectionally 
coupled in a master-slave configuration including appropriate delay terms, the slave can 
predict the trajectory of the master [2j. Remarkably, it has been proved that this kind of 
synchronization is stable and robust even in the presence of an external forcing acting upon 
both systems. Besides the many theoretical papers [S'-'S], anticipated synchronization has 
been experimentally observed in electronic circuits [U HDl E] and laser systems [121 US] and 
it has been proposed as a mechanism to control dynamical systems [HHTH] , competing with 
more traditional techniques as the Ott, Grebogi and York [T7j or the Pyragas [18] meth- 
ods. The idea behind the control technique using anticipated synchronization, named as 
the predict-prevent control method [15] , lies in the use of the information obtained from the 
slave output to prevent unwanted behaviors in the master system. 

In this paper we perform an extensive numerical study of the predict-prevent control 
method for excitability using FitzHugh-Nagumo model systems [T9l 120] as prototypic exam- 
ples. We show that it is possible to control the master system by monitoring the output of 
the slave by using this method. We study the robustness of the control for different types 
of external forcing as well as with respect to the parameter mismatch between master and 
slave. 

Although anticipated synchronization has been described before for excitable systems [Sj 
lll|211[22], a detailed characterization of the anticipated synchronization region in parameter 
space and its dependence with the forcing is still lacking, specially when using different 
forcing terms in the master and the slave. Therefore, we devote section [IT] to the description 
of the model equations as well as the characterization of the anticipated synchronization 



regions. In section III we show how to control the master system by using the output of the 



slave system. In section IV we analyze the influence of parameter mismatches. Finally, in 
section M we present the summary and conclusions. 
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II. EQUATIONS AND FORCING SCHEMES 



We consider two FitzHugh-Nagumo systems in the presence of external forcing terms and 
unidirectionally coupled in a master-slave configuration. The (dimensionless) equations are 

m 

xi{t) = - a)(xi - 1) - X2 + (1) 

X2(t) = e{xi - 6x2), (2) 

yi{t) = -yiivi - a){y^ - 1) - 7/2 + hit) + K[xi{t) -yiit- r)], (3) 

mit) = e{yi~by2), (4) 

where {xi,X2) are the dynamical variables of the master system and {yi,y2) are the cor- 
responding ones for the slave system; a = 0.139, b = 2.54 and e = 0.008 are constant 
parameters (as e ^ 1 there is a separation of time scales between the fast, xi,yi, and the 
slow, X2,y2, variables). and hit) represent external forcings. In Eq. ^ the term 

K,Xi{t) represents the unidirectional coupling from the master to the slave, while —nyi{t — T) 
is a feedback term delayed in time by an amount r (k controls the strength of both the 
coupling and the feedback). Note that variables Xi and yi are coupled in such a way that 
the dynamics of the master influences, but is independent of, the dynamics of the slave. 

Let us first describe briefly the main characteristics of the dynamics of the master system. 
Without external forcing, there exists a fixed point at the origin xi = X2 = 0. This fixed 
point is stable against small perturbations induced by the external forcing /i(t), but if 
the perturbations exceed a threshold value, then the system returns to the fixed point 
by a large excursion in phase space, a so-called spike or pulse. If, on the other hand, 
h{t) = Iq, a small constant, the stable fixed point is slightly different from zero. More 
specifically, for Jq < 0.035 (for the set of parameters a, b, e noted before) the fixed point is 
stable (and excitable behavior can occur), whereas for Jq > 0.035 the dynamics is oscillatory. 
Of course, the case of a constant forcing is not very interesting and we will consider that the 
forcing acting on the master system can be decomposed as hit) = Iq + ^i(t), being ^i(t) a 
random function of time and we take throughout the paper Jq = 0.03, below the excitability 
threshold. We have considered two possibilities: (1) ^i{t) is white noise with zero mean 
and correlations (^i(t)^i(t')) = D6{t — t'); (2) ^i(t) = Dj^k^it — tk) is a sum of impulses 
at random times tk such that the time differences tk+i — tk are distributed according to an 
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exponential distribution of mean value A. In both cases we call D the noise intensity. As a 
consequence of the random forcing, the system displays pulses at random times, see Fig. [1} 

Let us now consider the dynamics of the slave system. For constant common forcing 
{hit) = /2(t) = /o), there is a solution of the previous equations in which yi(t) = xi(t + r) 
and 2/2 (^) = X2{t + T). This remarkable solution, first found by Voss [2J, shows that the slave 
anticipates (i.e. predicts) by an amount of time r the dynamics of the master. Our intention 
is to use this anticipation property of the slave to influence the dynamics of the master in 
order to suppress all unwanted pulses. 

For non-constant Ii{t) it is no longer true that the anticipated manifold is an exact 
solution even in the case of common forcing. We will consider that the forcing on the slave 
is either constant, l2(t) = Iq, or it can be decomposed as l2(t) = /q + ^2(t) with C,2(t) a 
white noise of intensity D. In some cases we will take ^2(^) = ^i{t) and in others that C,2{t) 
and ^i(t) are independent random processes. We have shown in previous work |3] that, in 
this more general forcing scenario, it can exist nevertheless a region in the parameter space 
{k,t) such that it holds that yi{t) ~ + r) and, in particular, that the slave can display 
pulses that anticipate by a time approximately equal to r the pulses of the master. In the 
following, we analyze the details of this anticipated synchronization region in four cases: 

(i) identical random white noises ^i(t) = C,2{t), 

(ii) white noise ^i(t) in the master and constant forcing Iq in the slave, 

(iii) independent white noises ^i(t) and ^2(t), and 

(iv) a sum of impulses for C,i{t) in the master and a constant forcing Iq in the slave. 

In each case, we fix the values of Iq and D and study the region of parameter space (k, r) in 
which the slave anticipates correctly the pulses of the master. In Fig. [T]we show an example 
of a trajectory and a detail of the anticipated synchronization in Fig. [2] 

The region in parameter space in which anticipation of pulses is possible is shown in Fig. 
[3] for the different forcings (i)-(ii)-(iii)-(iv) described before. In each case we have set fixed 
values of Jo, D and quantified the quality of the anticipated synchronization by the following 
measure: pi = \n{xi) —nc{yi)\/n{xi) being n{xi) the number of pulses in the master system 
observed in a given (large) time interval and nc{yi) the number of those pulses which are 
correctly predicted by the slave system. We have considered that a pulse at a time t in the 
master is correctly predicted if there is another pulse in the slave occurring at a time t' {t 
(resp. t') are taken when the master (resp. slave) crosses the value 0.6) such that At = t — t' 
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satisfies < At < 20 (this time being about half the width of a typical pulse). As it turns 
out that the slave might display some extra, spurious, pulses with no correspondence in the 
master, it is also necessary to introduce p2 = \n{yi) —n{xi)\/n{xi) , being n{yi) the number of 
pulses observed in the same time interval in the slave system. Furthermore, as we are looking 
for a real anticipation between the pulses, we need to know that the average anticipation 
time is close to r and, hence, we introduce the third measure ps = |(At) — r|/r. Perfect 
anticipated synchronization implies pi = P2 = Pa = 0. The anticipated synchronization 
region in Fig. |3] signals the values of (k, r) for which pi < 0.1, p2 < 0.1 and pa < 0.2. In 
agreement with the results of reference [IS] , in which case (i) was considered using a different 
measure of synchronization, we find that for too large or too small coupling k, anticipation is 
lost. Moreover, a too large delay time r also prevents anticipated synchronization of pulses 
to occur. In the next section, we will consider parameter values in which anticipation does 
occur and we will devise a method that will allow us to suppress in an efficient way the 
pulses in the master system. 



III. CONTROL 



As stated before, our goal is to suppress the appearance of the random pulses in the master 
system, consequence of an unavoidable random forcing Ii{t). To this end, we consider first a 
simple direct control procedure that consists in reducing the amplitude of the xi{t) variable 
whenever it surpasses a threshold value xq > 0. This simple scheme is such that if to is the 
time at which xi(to) = crosses the threshold value, then we set xi(t = tg) = xi{to) — p, 
being p the amplitude of the control 'kick'. This is equivalent to adding a term of the form 
— p5(t — to) to Eq. ([T|. In this direct control method, we have to be careful about the precise 
value of Xq. If Xq ~ then the condition x{t) = xq might not lead necessarily to a pulse 
and we would be applying a control condition when it is not needed. If, on the other hand, 
Xq is very large, the pulse will be already too developed and it will be necessary to apply a 
strong control amplitude p in order to suppress it. 

In the following, we consider a modification of this simple direct control method. Follow- 
ing the ideas developed in reference [T5], we show that the use of a predict-prevent control 
scheme, based upon the slave variable crossing the threshold value, yiit'o) = xq, is more effi- 
cient in the sense that it allows to suppress more pulses with a smaller control amplitude p. 
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Likewise, the use of a value of p as small as possible is important in order not to introduce an 
uncontrolled perturbation to the dynamics of the master, as large values of p might induce 
the appearance of additional pulses after the suppressed ones. 

In Fig. |4]we plot the resulting time series for the variable Xi{t) after applying the control 
procedure just described in case (ii): white noise applied to the master and a constant 
forcing to the slave. It can be clearly seen that, for the same threshold value xq and control 
amplitude p, a more efficient control procedure (i.e. a larger fraction of suppressed pulses) 
occurs when using the predict-prevent control scheme based upon the slave variable than 
the one based on the direct control method of the master. The reason is obvious, as the 
pulses of yiit) precede by a time approximately equal to r the corresponding pulses of xi{t), 
the control is applied earlier in time, when the pulse in the xi{t) is not so well developed 
yet and it is easier to suppress. In all cases, and in order to avoid spurious and repeated 
control kicks, we have set a recovery time tree = 100 (this is a value larger than the average 
time- width of a pulse), such that two consecutive correcting kicks can not be applied in a 
time shorter that tree- 
In Fig. [5] we quantify the error of the control scheme by plotting the fraction of pulses 
that were not successfully suppressed as a function of the control amplitude p. As mentioned 
above, for large p the control procedure might lead to the appearance of new pulses, but, in 
any event, we see that the number of remaining pulses is a decreasing function of p. This 
figure shows that a substantial improvement in the suppression of pulses is obtained when 
using the predict-prevent control method based on the slave system. This is the main result 
of this paper. 

To show the robustness of the control scheme based on the slave system, and in order 
to cover a wider range of possible experimental situations, we show in the same figure the 
fraction of not suppressed pulses for the whole set of forcing schemes (i)-(ii)-(iii)-(iv) defined 
before. In all much better reduction with respect to the direct control using only 

the master is achieved. 

We have also studied the effect of a time lag tR between the crossing of the signal control 
with xo and the application of the control at xi. As it can be seen in the inset of Fig. [5] for 
forcing case (ii), as tji increases, less pulses are suppressed, both in the direct and the predict- 
prevent control procedures. These results imply that for a successful pulse suppression the 
time between the crossing of the signal and the application of the control must be as short 
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as possible. We have observed similar results for the other forcing cases. 

IV. PARAMETER MISMATCH 

As it is very unlikely that one can produce a perfect copy of the master system, an 
important issue concerning the control based on a master-slave configuration is how robust 
it is upon differences in parameter mismatch between the two systems. In fact, an experiment 
using an electronic implementation of the FitzHugh-Nagumo neurons has been carried out 
in reference [I5] with the result that the control procedure can be carried out safely with 
real, non-identical, systems. 

We have studied the effect of small variations of the parameters a, b and e in the slave 
system, Eqs. ([3]|4]), in order to analyze how the control of the master dynamics is affected. 
We have first studied the effect of each one of these parameters separately by changing about 
5% its value. As a consequence the fraction of remaining pulses changes by: 20% in the case 
of varying a, 10% in the case of varying b, and 8% in the case of varying e. For values of a 
in the slave larger than in the master less pulses are controlled. A similar result applies to 
variations of e. However, the effect of b is different: for values of b larger in the slave than 
in the master, the fraction of remaining pulses decreases when using the predict-prevent 
control method, i.e, more pulses are controlled. The effect of a and b can be understood in 
terms of the stable fixed point: For different parameters in the slave and the master, the 
stable fixed points will be different. For larger values of a in the slave system, the fixed 
point moves towards the origin, so decreasing the excitability threshold. In contrast, for b 
larger in the slave the fixed point moves away from the origin, so increasing the excitability 
threshold and worsening the pulse control. 

In Fig. [6} the solid lines are the same than the ones shown in Fig. |5| where a, b and e 
in the slave equations are equal to those of the master equation. We have studied how the 
control dynamics is affected when changing at the same time the three parameters (a, b, e) 
in the slave equations. For up to a 5% (resp. 10%) of difference between the parameters of 
the master and the slave, the fraction of pulses that were not successfully controlled for xi 
are located between the two dashed (resp. dotted) lines. An specific example is to take the 
following parameters for the slave system: (a,6, e) = (0.1251,2.794,0.0072), i.e. a smaller, 
b larger and e smaller in the slave than in the master. In this case, less pulses remain, so 
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making the control procedure more effective (dotted line at the left of Fig. [6]). This shows 
that parameter mismatch does not necessarily worsen the results, as in some cases even 
a larger fraction of pulses can be successfully controlled. It is worth noting that we have 
observed similar results for the different forcing types that we have considered. 

V. CONCLUSIONS 

In this paper we have first characterized the anticipated synchronization region of two 
unidirectionally coupled Fitzhugh-Nagumo systems in a master-slave configuration for dif- 
ferent types of forcing terms, obtaining qualitatively similar regions for the different cases 
(i)-(ii)-(iii)-(iv) considered. This implies that anticipated synchronization is indeed a robust 
phenomenon even under different types of uncommon forcing for the master and the slave 
equations. 

Then we have performed a numerical study of the predict-prevent control method. For 
parameter values inside the anticipated synchronization region, we have applied two control 
schemes in order to suppress the pulses of the master: the direct control using the master 
output and the predict-prevent control using the slave output. We have obtained that the 
predict-prevent control is more efficient than the direct control using only the master. This 
statement is true for a variety of random forcing terms, both common and uncommon to 
the master and the slave. 

Finally, we have obtained that a perfect matching between the parameters of the master 
and the slave is not necessary for the control to be efficient and, in fact, a slight parameter 
mismatch can, in some cases, lead to a better control. 

The results obtained in this work are a clear indication of the robustness of the proposed 
predict-prevent control method and opens the door to more general experimental implemen- 
tations, in other physical and biological systems, than the ones carried out previously [15]. 
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FIG. 1: (Color online) Trajectories xi{t) and yi{t) coming from a numerical integration of Eqs. 
([T][4]). The forcing in the master is always Ii{t) = Iq + ^i{t) being ^i(t) white-noise, and the corre- 
sponding trajectory is plotted in the bottom panel (black line). We also plot the slave trajectories 
when different forcing schemes (i)-(ii)-(iii) (see the main text) are applied to the slave. Parameters: 
a = 0.139, b = 2.54, e = 0.008, k = 0.4, r = 2.4, Iq = 0.03, D = 2.45 • 10"^ 
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FIG. 3: (Color online) The dots show the anticipated synchronization region defined as the set of 
values of (k,t) for which pi < 0.1, p2 < 0.1 and ps < 0.2 in the four forcing cases (i)-(ii)-(iii)- 
(iv) explained in the main text. Other parameter values as in Fig. [l] except in (iv): Iq = 0.032, 
D = 0.05, A = 500. 
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FIG. 4: (Color online) Time traces of the master variable xi{t). In panels A, left, we show 
the pulse suppression after applying a direct control ('kick') of magnitude p whenever the master 
variable xi{t) crosses the threshold value xq = 0.3. As p increases the number of pulses successfully 
suppressed increases with respect to the no-control case, p = 0. In the right panels, B, we show 
that the reduction is more effective if we use a control scheme based upon the slave variable yi{t) 
crossing the same threshold value. The results come from a numerical integration of Eqs. (l][4) in 
the forcing scheme (ii): a random forcing Ii{t) = /q + '^i(0; being ^i(t) Gaussian noise of intensity 
D, in the master and a constant forcing /2(i) = -^o in the slave. Parameters are: a = 0.139, 
b = 2.54, e = 0.008, Iq = 0.03, D = 2.45 • 10-^ k = 0.4, r = 2.4 (k and r correspond to the black 
dot in Fig. |3j). The recovery time after which no other correction can be applied is tree = 100. 
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FIG. 5: (Color online) Fraction of remaining pulses of xi as a function of the control parameter 
p. The lines with + symbols correspond to the direct control scheme based only upon the master 
system for forcing schemes (i)-(ii)-(iii) (solid) and (iv) (dot-dashed). Same parameter values than 
in Figs. [3] and |4j The different lines at the left of the main graph correspond to the predict- 
prevent control using the slave in the different forcing schemes: (i), red dashed; (ii), blue solid; (iii) 
green dotted; and (iv), purple dot-dashed. In the inset, we show the effect of a reaction time tn 
that elapses between the control variable crosses the threshold value xq and the corrective impulse 
is applied for case (ii) and p = 0.04 using direct control on the master (solid line with -|-) or 
predict-prevent control on the slave (solid line). 
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FIG. 6: (Color online) In order to show the effect of parameter mismatch, we plot the fraction 
of remaining pulses of function of p in forcing case (ii) using the predict-prevent con- 

trol method. The solid line corresponds to identical parameters in master and slave: (a, b, e) = 
(0.139, 2.54, 0.008). The values between the dotted (resp. dashed) lines correspond to the parame- 
ters in the slave varying up to 10% (resp. 5%) with respect to those of the master. The solid line 
with + symbols is the same than in Fig. [5] and it corresponds to the direct control method and is 
plotted here for comparison. Other parameters as in Figs. [3] and [4} 
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